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Abstract. We determine the mean velocity dispersion of six Galactic outer halo globular clusters, AM 1, Eridanus, 
Pal 3, Pal 4, Pal 15, and Arp 2 in the weak acceleration regime to test classical vs. modified Newtonian dynamics 
(MOND). Owing to the nonlinearity of MOND's Poisson equation, beyond tidal effects, the internal dynamics of 
clusters is affected by the external field in which they are immersed. For the studied clusters, particle accelerations 
are much lower than the critical acceleration ao of MOND, but the motion of stars is neither dominated by internal 
accelerations (a< S> a e ) nor external accelerations (a e 3> ai). We use the N-body code N-MODY in our analysis, 
which is a particle-mesh-based code with a numerical MOND potential solver developed by Ciotti, Londrillo, 
and Nipoti (2006) to derive the line-of-sight velocity dispersion by adding the external field effect. We show that 
Newtonian dynamics predicts a low- velocity dispersion for each cluster, while in modified Newtonian dynamics the 
velocity dispersion is much higher. We calculate the minimum number of measured stars necessary to distinguish 
between Newtonian gravity and MOND with the Kolmogorov-Smirnov test. We also show that for most clusters 
it is necessary to measure the velocities of between 30 to 80 stars to distinguish between both cases. Therefore the 
observational measurement of the line-of-sight velocity dispersion of these clusters will provide a test for MOND. 
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1. Introduction 

Observable matter in galaxies and in clusters of galax- 
ies cannot produce sufficient gravity to explain their dy- 
namics. Cold dark matter (CDM) scenarios or alterna- 
tive theories of gravitation are therefore invoked to re- 
solve the problem. Nowadays, the CDM hypothesis is the 
dominant paradigm. This hypothetical matter does not 
interact with electromagnetic radiation and only shows 
its presence through its gravitational interaction. Even 
though the dark matter hypothesis has successfully ex- 
plained the internal dynamics of galaxy clusters, grav- 
itational lensing, and the standard model of cosmol- 
ogy within the framework of general relativity (GR) 
(Spergel 2003), much experimental effort has failed to 
yield a detection of dark matter particles. Moreover the 
results of high-resolution simulations of structure forma- 
tion do not reproduce some observations on galactic scales, 
such as the central structures of rotation curves, the preva- 
lence of low bulge-to-disc ratios, and the numbers and 
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spatial distribution of the subhalos (|Klypin et al. 1 999 
IMoore et al. 19991 IMetz et al. 20081 |Kroupa et al. 2010p . 
Even the ability of the dark matter theory to account for 
the Tully-Fisher and Freeman relations is controversial 
(|Bosch fc Dalcanton 20001 IGovernato et al. 2010|) . These 
shortcomings have not led to the rejection of the theory 
only because on galactic scales baryons are at least non- 
negligible contributors to the mass density, consequently 
simulations that include the complex physics of star for- 
mation are essential for reliable predictions. Currently 
such simulations are still at an experimental stage and 
are usually substituted by "semi-analytic" arguments that 
have weak theoretical underpinnings. 

One of the alternative theories to CDM is the so- 
called modified Newtonian dynamics (MOND) theory, 
which was originally proposed by Milgrom (1983) to ex- 
plain the flat rotation curves of spiral galaxies at large 
distances by a modification of Newton's second law of 
acceleration below a characteristic scale of ao ~ 1.2 x 
10 -10 ms -2 = 3.6 pc/(Myr) 2 without invoking dark mat- 
ter ( [Bekenstein fe Milgrom 1984[ ). In MOND, Newton's 
second law is modified to /z(a/do)a = ajy + V x H, where 
p is the mass-density distribution, is the Newtonian 
acceleration vector, a is the MONDian acceleration vec- 
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tor, a — |a| is the absolute value of MONDian accel- 
eration, (i is an interpolating function for the transition 
from the Newtonian to the MONDian regime, which runs 
smoothly from = i at i < 1 to n(x) = 1 at 

x 1 (Bekcnste in fe Milgrom 1984 ). Different interpo- 
lating functions have been suggested, such as the simple 
function, [i{x) = x/(l + x) QFamaey fe Binney 20 05) and 
the standard interpolation function, = x/y/l + x 2 

(Milgrom 1983). Because the simple function fits galac- 
tic rotation curves better than the standard function 
(Gentile, Famacy & de Blok 2010 1, in this paper we use 
the simple function. The value of the curl field H de- 
pends on the boundary conditions and the spatial mass 
distribution and vanishes only for some special symme- 



tries (Bekenstein fc Milgrom 1984). The non-linearity of 



the MOND field equation leads to difficulties for standard 
N-body codes and makes the use of the usual Newtonian 
N-body simulation codes impossible in the MOND regime. 

It has been shown that on galac- 
tic scales MOND can explain many phe- 
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alized to a general-relativistic version (Bekcnstein 2004, 
ISanders 20051 |Zlosnik, Ferreira fc Starkman 2007[ 

Mil grom 2009[ ), making it possible to test its predic- 
tions for gravitational lensing. Dynamics of galaxies in 



clusters ( |Sanders fc M cGaugh 2002) and the merging 
of galaxy clusters, where the baryonic mass is clearly 
separated from the gravitational mass ([Clowe et al. 2006)) 
cannot be completely explained by MOND without 
invoking some kind of hot dark matter, perhaps in the 
form of massive (active or sterile, 2 to 11 eV) neutrinos 
QAngus et al. 2006||Xngus et al. 2010 1 . 

In order to decide whether MOND is a comprehensive 
theory to explain the dynamics of the universe, it is desir- 
able to study MOND for objects in which no dark matter 
is supposed to exist and where the characteristic acceler- 
ation of the stars is less than the MOND critical acceler- 
ation parameter ao- Globular clusters (GCs) are a perfect 
candidate since they are the largest virialized structure 
that do not contain dark matter (Moore 1996). 

In the distant halo of our Milky Way there exist several 
low-mass GCs where both internal and external accelera- 
tions of stars are significantly below the critical accelera- 
tion parameter ao of MOND. Because GCs are assumed 
to be dark-matter-free, if MOND is true, the motions of 
stars must deviate from the standard Newtonian dynam- 
ics. It has been proposed by Baumgardt (2005) that some 
of these distant Galactic GCs are perfect tools to test grav- 
itational theories in the regime of very weak accelerations. 
For MOND, the internal velocity dispersion among the 
stars in these clusters would be significantly higher than 
in Newtonian dynamics. 



The mean velocity dispersion of stellar systems for 
the two extreme cases of internal (a^ 3> a e ) or external 
(a e 3> a<i) field domination have been derived analyti- 



cally by Milgrom dMilgrom 19861 Pil grom 1994| ), assum- 
ing that the systems are everywhere in the deep-MOND 
regime (a e ,a; <C ao)- Many systems that can be used to 
test MOND are not completely either internally or exter- 
nally dominated. Globular clusters or dwarf galaxies of the 
Milky Way for example have internal and external acceler- 



ations that are of the same order ( Baumgardt et al. 2005 ), 
consequently one has to determine the velocity dispersion 
numerically for intermediate cases. Sollima and Nipoti 
(2009) constructed self-consistent, spherical models for 
stellar systems in MOND, neglecting the external field 
effect and presented a dynamical model for six galactic 
globular clusters. The presence of the external field ef- 
fect breaks the spherical symmetry and validity of their 
model. Recently Haghi et al. (2009, hereafter HBK09) in- 
vestigated the dynamics of star clusters by numerically 
modeling them in MOND, assuming circular orbits. They 
performed N-body simulations and presented analytical 
formulae for the velocity dispersion of stellar systems in 
the intermediate MOND regime, which are useful for a 
comparison with observational data of several GCs and 
dSph galaxies (for details on the numerical calculations 
see HBK09). In a follow-up paper, Jordi et al. (2009) 
determined the velocity dispersion (using 17 stars) and 
mass-function slope of Pal 14 and showed that MOND 
can hardly explain the low- velocity dispersion of this sys- 
tem. However, Gentile et al. (2010) showed that with 
the currently available data, the Kolmogorov-Smirnov 
(KS) test is still unable to exclude MOND with a suffi- 
ciently high confidence level. Moreover, the low density 
of Pal 14 suggests that binary stars may be an impor- 
tant issue for interpreting its measured velocity disper- 
sion QKupper fc Kroupa 2010[ ), and the true velocity dis- 
persion of Pal 14 could be much lower than the value re- 
ported by Jordi et al. (2009), thereby possibly posing an 
even larger challenge for MOND, but also for Newtonian 
gravity and for any understanding of the dynamics of this 
object as being in equilibrium. 

In this paper we calculate the prediction of MOND 
and Newtonian dynamics on the velocity dispersion of six 
other distant clusters of the MW (Table 1). In order to 
see the pure MONDian effects, we concentrate on sys- 
tems in which the tidal radius is much larger than the 
gravitational radius Q and therefore tidal effects are unim- 
portant. In other words, this paper provides the basis for 
further observational efforts. The measurements of a low- 
Newtonian) velocity dispersion would mean that MOND 
in its present form is in severe trouble and that globular 



1 The gravitational radius is a measure of the size of the 
system and is related to the mass and potential energy as given 
in Eq. 2-132 in Binney and Tremaine (1987). In many stellar 
systems the gravitational radius can be approximated by the 
three-dimensional half-mass radius as r 9 = 1.25rh if the 
assumption of virial equilibrium is valid. 
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clusters do not possess dark matter. In contrast, a high- 
velocity dispersion would either favor MOND or could be 
a hint to the existence of dark matter in globular clusters 



dBaumgardt et al. 2009 1. 

The paper is organized as follows: In Section 2 we give 
a brief review of the external field effect (EFE) in MOND. 
The simulation setup is explained in Section 3. The nu- 
merical results for six clusters are discussed in Section 4. 
We present our conclusions in Section 5. 

2. External field effect in MOND 

In classical Newtonian dynamics, a uniform external field 
does not affect the internal dynamics of a stellar system. 
In other words, in the frame of the system relative motions 
of objects are the same as in an isolated system. 

In MOND, the situation is entirely different. 
Owing to the non-linearity of Poisson's equation, 
the strong equivalence principle (SEP) is violated 
( |Bekenstein &: Milgrom 1984[ ), and consequently the in- 
ternal properties and the morphology of a stellar sys- 
tem are affected by both the internal and external field. 
This so-called external field effect (EFE) significantly af- 
fects non-isolated systems and can provide a strict test for 
MOND. 

The EFE is indeed a phenomenological requirement 
of MOND and was postulated by Milgrom (1983) to ex- 
plain the dynamical properties of open clusters in the 
MW that do not show MONDian effects. The EFE al- 
lows high-velocity stars to escape from the potential of 
the Milky Way ( |Famaey et al. 2007| IWu et al. 2007| and 
implies that rotation curves of spiral galaxies should fall 
where the internal acceleration becomes equal to the ex- 
ternal acceleration (jGentile et al. 20071 IWu et al. 2008[) . 

For a star cluster with a density distribution p c , which 
is embedded in a host galaxy with a density distribution 
p ex ti the acceleration of stars in the cluster satisfies the 
modified Poisson equation 



V.[/i( V$ = 4nG( Pc + p ext ), 

a 



(1) 



where V$ is the MONDian potential generated by the 
total matter density. For star clusters or dwarf galaxies 
far out in the halo of the Milky Way the local density 
of the Milky Way is negligible ( i.e. p ext « p c ). One 
way to solve Eq. 1 is then to assume that the total ac- 
celeration is the sum of the internal a.; and the external 
a e acceleration, which both satisfy the modified Poisson 
equation as (|Bekenstein fe Milgrom 1984|IWu et al. 20071 
IWu et al. 20081 |Angus 2008| ) 



V.[ M ( 



|a e + aj| 

a 



)(a e +a ; )] =AttG Pc 



(2) 



In Eq. 2, the direction of the external acceleration is 
important for the dynamics of stars in the clusters. For ex- 
ample, for two stars the acceleration would be different if 
a e is parallel or anti-parallel to their internal acceleration. 



Several attempts have been made considering the EFE. 
Wu et al. (2007) simulated an isolated object with a 
static MONDian potential solver, but changed the bound- 
ary condition on the outermost grid point to be nonzero. 
Famaey et al. (2007) estimated the escape velocity of 
galaxies in MOND and assumed /x(|a e + ai|/ao)ai = a^. 
As a first order test they replaced la^+ail with (a e +cti), or 
yal + of . This is an approximation, since they neglected 
the direction of a e (i.e. a possible angular difference be- 
tween ai and a e ). 

In order to study the evolution of the systems we need 
to calculate the acceleration at each step by adding the 
constant external gravitational field to the internal accel- 
eration inside the /i-function, considering the different an- 
gle between external and internal acceleration for all stars 
throughout the evolution. 

3. Simulation setup 

We performed a large set of N-body simulations of star 
clusters with the N-MODY code, which has been devel- 
oped by Ciotti et al. (2006). N-MODY is a parallel, three- 
dimensional particle-mesh code for the time-integration of 
collision- less N-body systems ( ]Londrillo fc Nipoti 2 009). 
The code numerically solves the non-linear MOND field 
equations, which can be used to perform numerical exper- 
iments in either MONDian or Newtonian dynamics. The 
potential solver of N-MODY is based on a grid in spher- 
ical coordinates and is best suited for modeling isolated 
systems. N-MODY uses the leap-frog method to advance 
the particles. The code and the potential solver have been 
presented and tested by Ciotti et al. (2006) and Nipoti et 
al. (2007). 

In the present study we used a spherical grid (r, 0, ip) 
made of N r x Ng x N v = 128 x 128 x 128 grid cells for the 
integration. The total number of particles is N p = 10 5 . 
The details of the scaling of the numerical MOND models 
and code units are discussed in Nipoti, Londrillo & Ciotti 
(2007). 

Our treatment of the external field effect is based on 
the idea that the star clusters we study are much smaller 
than they are distant from the Milky Way, so the external 
field is nearly constant over the cluster area. We carry 
out the simulations in a rotating reference frame centered 
on the cluster. We also assume that the clusters are on 
a circular orbit, so the external field of the Milky Way is 
also constant with time in the clusters frame. 

We assume that tidal forces arising from a gradient 
of the external field and the Coriolis acceleration caused 
by the rotating reference frame can be neglected. This is 
justified as long as the size of the cluster is much smaller 
than the tidal radius. We then solve Eq. 2 numerically 
with N-MODY, using a constant acceleration a e as the 
boundary condition. 

In order to include the EFE for non-isolated systems, 
we applied several changes to the N-MODY code. The 
changes were encoded in the source file mond-lib.f90, 
which contains all relevant routines implementing the 
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Cluster name 


Rh [pc] 


R G [kpc] 


a e [ao] 


M C [10 3 M Q ] 


<jm [kms -1 ] 


<7jv[kms 1 ] 


-^V min 


AM 1 


24 


123.2 


0.086 


12.6 


1.50 


0.56 


25 


Eridanus 


14.2 


95.2 


0.113 


18.6 


1.80 


0.90 


40 


Pal 3 


24 


95.9 


0.112 


31.6 


1.97 


0.87 


30 


Pal 4 


23.2 


111.8 


0.096 


42.6 


2.30 


1.10 


35 


Pal 15 


21.2 


37.9 


0.283 


26.3 


1.42 


0.88 


80 


Arp 2 


21.5 


21.4 


0.504 


21.8 


1.07 


0.80 


150 



Table 1. Globular clusters modeled in this paper. The half- mass radius and galactocentric distances, Rq, are taken 
from Harris (1996). The values of the external acceleration are calculated from a e = -^/GMclq/Rg with M = 1.2 x 
10 n MQ for all galactocentric distances. Cluster masses, M c , were calculated from the absolute V-band luminosities 
by assuming a stellar mass-to-light ratio of MjLy = 2, which is close to the measured average mass-to-light ratio of 
galactic globular clusters (Mieske et al. 2008). um and <tn are the corresponding MONDian and Newtonian values of 
the velocity dispersion, respectively. The last column is the minimum number of stars necessary to obtain P < 0.05. 



MOND potential solver. The Mondian potentials and ac- 
celerations are assigned to stars in subroutines: " gmond" 
and " mond-acc" , and we placed the vector of a e in the in- 
terpolating function for all stars, i.e., a e is added to a, via 
vector summation and the code solves Eq. 2 by iteration. 
At each step of potential solving, we added the constant 
external field for all grid points according to Eq. 2. This 
is a first approximation that allows us to focus on the ef- 
fects of a constant external field, which in turn allows us 
to deal with the intermediate MOND regime. We would 
like to stress that our numerical solution agrees very well 
with analytic ones for the extreme cases <C a e <C ao, 
a e <C a,i <S ao and a e ,a,i ^> ao (HBK09). 

Our method is a first attempt to take into account 
the constant EFE in N-MODY, in order to avoid solving 
the modified Poisson equation for the large area including 
both galaxy and cluster, which is impossible with the cur- 
rent version of N-MODY. This can only be studied with a 
realistic, high-resolution MOND simulation that includes 
the density distribution causing the external field. 

In this work we start from Newtonian equilibrium 
Plummer models and produce MONDian equilibrium ini- 
tial systems. In order to have a MONDian equilibrium 
initial system, we increased the initial velocity of the par- 
ticles. This method is useful when the external field is 
important. The details of the method are described in 
HBK09. 



4. Modeling the distant globular clusters 

We calculated the velocity dispersion for the Galactic 
globular clusters mentioned in Table 1. These clusters are 
best suited for testing MOND. A large set of dissipation- 
less N-MODY computations with MONDian equilibrium 
initial conditions were performed for stellar systems that 
are embedded in the outer MW halo and affected by dif- 
ferent values of the external field. Clusters in Table 1 are 
generally far out in the Galactic halo so that the external 
acceleration caused by the MW is small. In order to pro- 
duce different internal acceleration regimes, we changed 
the cluster mass from 500M Q to 10 7 M Q and assume the 
half-mass radius to be constant. The models are evolved 



for several crossing times to reach the equilibrium state. 
The EFE on the predicted line-of-sight velocity dispersion 
for non-isolated stellar systems with different internal ac- 
celerations have been carried out in HBK09. 

For the GCs listed in Table 1, observational efforts 
are underway to determine their velocity dispersion and 
to constrain their mass (Grebel et al. 2009). Because the 
half-mass relaxation times of most clusters in Table 1 are 
on the order of a Hubble time or even larger, dynamical 
evolution does not play an important role for these clus- 
ters, therefore we can assume that the half-light radius 
is equal to the half-mass radius. The projected half-mass 
radii, Rhp, and galactocentric distances, Rq, are taken 
from Harris (1996). The projected half-mass radii relate 
to the three-dimensional half-mass radii as Rh P = 7-R/i 
with 7 w 0.74. For all clusters in Table 1, the tidal radii 
are larger than the half-mass radii by a factor of 5 to 10, 
which means that tidal effects play no significant role for 
the internal dynamics of the clusters in our sample. We 
perform numerical modeling to obtain the mean velocity 
dispersion. 

Because the cluster masses are not known from ob- 
servations, it is useful to calculate line-of-sight velocity 
dispersions for different values of cluster mass. If the ve- 
locity dispersion is determined observationally, one can 
constrain the cluster mass. In Fig. 1 the resulting global 
line-of-sight velocity dispersion as a function of mass is 
plotted for the six stellar clusters and is compared with the 
Newtonian results. The black line shows the Newtonian 
prediction for the velocity dispersion, which we calcu- 
late using a L os,N = O.Z6y/GM/Rh (|Haghi et al. 2009[). 



The red line shows the numerical calculation of the ve- 
locity dispersion in the MONDian regime. For each clus- 
ter the external acceleration of the Galaxy is given in 
Table 1. All clusters are in the intermediate regime, for 
which there is no analytical prediction in MOND. In the 
low-acceleration region, the numerical solutions show a 
considerable relative difference with Newtonian results. 
As the mass increases, the internal acceleration grows 
and then gets close to the Newtonian results (at about 
M — 10 6 M Q ). The vertical dashed lines show the ex- 
pected mass for each cluster calculated from the absolute 
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3.5 4 4.5 



Eridanus 

3.5 4 4.5 




3 3.5 4 4.5 
log (Mass) [sollar mass] 



Fig. 1. Line-of-sight velocity dispersion for clusters listed in Table 1 for various masses as found by N-MODY. In 
order to compare with the real cluster, the half-mass radii of all models are fixed. For low masses, which mean low 
internal accelerations, the prediction of Newtonian dynamics differs from the numerical solution. As the mass increases 
(at about 10 6 Mq), the internal acceleration grows and the system enters the Newtonian regime and the numerical 
solutions get close to the Newtonian prediction. In order to focus on a reasonable range of mass for these clusters, 
we show the results in the range of 500 — 2 x 10 5 . Dashed vertical lines show the expected cluster mass, calculated 
from the absolute luminosity by assuming a stellar mass-to-light ratio of M/Ly = 2. This value agrees with observed 
mass-to-light ratios of galactic globular clusters (Mieske et al. 2008) 



luminosities by assuming a stellar mass-to-light ratio of 
MIL = 2 (jMieske et al. 2008|) . It should be noted that in 
the simulations we do not assume any M/L, but we vary 
the mass and therefore the M/L ratio of the clusters. The 
velocity dispersions corresponding to this mass for both 
cases are indicated in Table 1. The M/L ratio of 2 is only 
used in Table 1 to illustrate there is a MOND effect, which 
makes a difference for the studied clusters. Our predictions 
for the MONDian case are lower than those by Baumgardt 
et al. (2005). The absolute differences between MONDian 
and Newtonian results are largest for the clusters Pal 4, 
Pal 3, AM 1, and Eridanus. Therefore these clusters would 
be the best cases to test MOND. The same qualitative re- 
sults were also obtained by Sollima and Nipoti (2009), 
who have studied some of the clusters of our list (i.e., 
Eridanus, AMI, Pal 3, Pal 4, and Pal 14, which we have 
studied in HBK09). They have obtained the cluster's ve- 
locity dispersion profile. Their predictions for the velocity 
dispersion are generally higher than ours. For AMI for ex- 
ample, the mean value of the velocity dispersion is ~1.7 
kms -1 since the predicted MOND velocity dispersion pro- 
file is flat throughout the cluster area (Fig. 5 of Sollima 



and Nipoti (2009)). This is higher than the value of 1.5 
kms -1 we find, which is to be expected because they did 
not include the EFE. If we take the EFE into account, 
the systems tend to be in the quasi-Newotonian regime, 
where the gravity is weaker than in the deep-Mondian 
regime (HBK09). 

In MOND the external field introduces an anisotropy 
so that clusters could in principle become elongated. In 
order to test for a possible elongation of the velocity dis- 
persion in the direction of the external field (cr x ), we define 
the parameter a = 1 — CT^/(cr^ + al). Figure [2] shows that 
the line-of-sight velocity dispersion in the direction of the 
external field differs only randomly compared to that per- 
pendicular to the external field. This small a difference 
would be unobservable. 

In order to show the situation of the isotropy in the 
velocity dispersion, the anisotropy profile, /3(r) = 1 — 
(of /2of ) for one cluster is plotted in Fig. [5J where, a t 
is the tangential velocity dispersion and a r is the radial 
velocity dispersion. 

Figure [3] shows the distribution of radial velocities of 
the N fota ; = 10 5 stars in the velocity range -10 to +10 
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50 75 100 
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25 



_i_ 



50 75 100 



Fig. 2. Top panel: The ratio of the velocity dispersion in 
the direction of the external field compared to the perpen- 
dicular direction for Pal 3 as an example plotted against 
the radius. The system is isotropic throughout the cluster. 
Bottom panel: The anisotropy profile, j3, for Pal 3 as an ex- 
ample plotted against the radius. The system is isotropic 
throughout the cluster. 



km/sec for both Newtonian (red triangles) and Mondian 
(black dots) cases together with the best Gaussian fit to 
this distribution using the least square method. We cre- 
ated the Newtonian model with a mass of 177000M Q and 
the Mondian one with a mass of 32000M Q , so that the ve- 
locity dispersions are the same. The bottom panel shows 
the relative deviation, (N Newtonian-^ M and) /N total, where 
N is the number of stars in a velocity bin. Evidently the 
observed radial velocity distribution is well approximated 



0.1 



0.08 - 



» 0.06 h 

c 
c 
o 

'f 0.04 
0.02 h 



• Mond model 
A Newtonian model , 
Gaussian fit 



-10 



-5 



10 



0.02 



Z 0.01 - 



1-0.01 - 




Velocity (km/sec) 



Fig. 3. Upper panel: Radial velocity distribution of stars 
for a Newtonian (red triangles) and MONDian (black 
dots) system in which both have the same internal ve- 
locity dispersion. The smoothed line shows a Gaussian 
fitted to the distributions. Both distributions follow the 
Gaussian fit very closely. Lower panel: The relative differ- 
ence between Newtonian and Mondian velocity distribu- 
tions. There is a difference between the two models, but it 
is too small (less than 1 %) to be detected observationally. 



by this Gaussian. The Newtonian one is slightly different 
from the Mondian distribution, but the relative difference 
is very small (less than 1%), therefore it is impossible to 
find a detectable deviation observationally. 

5. Distinguishing between MOND and Newtonian 
gravity models 

Usually, it is not feasible to measure the velocity disper- 
sion of a large sample of stars in distant and low-mass 
globular clusters. If only a few stars are measured, the ob- 
served velocity dispersion has a large error, which makes 
it difficult to distinguish between MOND and Newtonian 
dynamics. Recently, Gentile et al. (2010) have shown that 
the small sample size for Pal 14 (i.e., radial velocity of 17 
stars were used to determine the observed velocity disper- 
sion) can only reject MOND with a low confidence level. 
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In order to see how many stars are necessary to reject 
MOND with high statistical confidence if the underlying 
velocity dispersion is Newtonian, we created artificial data 
sets of increasing sample size from the Newtonian simu- 
lations for each of the six clusters. For each sample size 
we produced 20 random realizations of a mock sample, 
and each time applied a KS test with the null hypothesis 
being that MOND is correct (thus comparing the mock 
data set from the Newtonian simulation with the veloc- 
ity dispersion of the MONDian simulation, ctmond)- We 
approximate the cluster mass from the V-band luminos- 
ity by assuming a stellar mass-to-light ratio of M/L = 2 
(vertical dashed line in Fig. 1). Figure 0] shows the mean 
P- value from the KS test as a function of the number of 
stars in the sample. The number of stars changes in the 
range of [10-160]. For any given cluster and a fixed num- 
ber of stars a 1-sigma range of P-values are shown as error 
bars. The minimum number of stars, iV m i„, which is nec- 
essary to achieve P < 0.05 (i.e. exclusion confidence of 
95%), are obtained for each cluster and are given in Table 
1. Measurement of N m i n stars should therefore be enough 
to rule out MOND if the cluster mass is simultaneously de- 
termined by independent measurements, possibly through 
e.g. star counts. 



6. Conclusions 

We compared global line-of-sight velocity dispersions of six 
distant low-density globular clusters of the Milky Way in 
MONDian and Newtonian dynamics and showed that they 
have a significantly higher velocity dispersion in MOND 
than the prediction of Newtonian dynamics. 

Using the N-MODY code, we obtained a large set of 
dissipationless numerical solutions for globular clusters 
with the MONDian initial conditions as end-products of 
the N-body computations by adding the external field. In 
order to produce different internal acceleration regimes, 
we changed the mass of the system over a wide range for 
each cluster. In addition, our results show that the clusters 
Pal 4, Pal 3, AM 1, and Eridanus are the best cases to test 
MOND owing to their larger absolute difference between 
MONDian and Newtonian velocity dispersions. These re- 
sults will allow us to test MOND more rigorously than was 
possible so far, and this will enable us to compare MOND 
with observational data in the future. 

Recently Sollima and Nipoti (2009) have performed a 
test for MOND by constructing self-consistent dynamical 
models for outer galactic clusters. These authors neglected 
the external field effect of the Milky Way, and hence ob- 
tained higher estimates for the velocity dispersion of clus- 
ters compared with our results. This difference is reason- 
able. Indeed, the external field in MOND is leading the 
system to move from the deep-MONDian regime toward 
the quasi-Newtonian regime, and thus to reduce the veloc- 
ity dispersion. In other words, the larger the external field, 
the smaller the internal acceleration, which implies a lower 



velocity dispersion for the cluster (Haghi et al. 2009). 



It should be noted that the clusters could be on ec- 
centric orbits, which means that the MOND predictions 
would be different because of the variation of the external 
field (in direction and amplitude) along the orbit. 

Using a KS test, we calculated the minimum number 
of stars that are sufficient to exclude MOND (under the 
hypothesis that these globular clusters are on circular or- 
bits) at the 95% confidence level. We found that between 
30 to 80 stars are necessary for most clusters to distin- 
guish between both cases. This number of stars can be 
observed with current 8m class telescopes. Additional ob- 
servational efforts to determine the velocity dispersions of 
these clusters and constraining the mass of the clusters 
by star counts would be highly important and provide a 
strict test of MOND. On the other hand, if MOND is 
the correct theory, these observations could be used to 
constrain the external field and consequently to put con- 
straints on the potential in which the systems are embed- 
ded. According to the anisotropy profile, the simulated 
systems are isotropic throughout. 
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Fig. 4. Mean P-value vs. number of stars in the sample for a number of random realizations of distributions of stars 
following a Gaussian with a standard deviation of cr/v for clusters listed in Table 1. The minimum number of stars 
necessary to reject MOND is defined where the mean P-value drops below 0.05. 
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